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We present a measurement of D° — D° mixing and CP violation using the ratio of lifetimes 
simultaneously extracted from a sample of D mesons produced through the flavor-tagged process 
D* + — > D°7r + , where D° decays to K^ir^, K~ K + , or 7r~7r + , along with the untagged decays 
D° -» T^tt* and D° K~ K + . The lifetimes of the CP-even, Cabibbo-suppressed modes K~ K + 
and 7r~7r + arc compared to that of the CP-mixed mode K^ir in order to measure yep and AY. 
We obtain y C p = [0.72 ± 0.18(stat) ± 0.12(syst)]% and AY = [0.09 ± 0.26(stat) ± 0.06(syst)]%, 
where AY constrains possible CP violation. The yep result excludes the null mixing hypothesis at 
3.3<r significance. This analysis is based on an integrated luminosity of 468 fb -1 collected with the 
BABAR detector at the PEP-II asymmetric-energy e + e~ collider. 



PACS numbers: 13. 25. Ft, 12.15.Ff, 11.30.Er 

I. INTRODUCTION 

Several measurements [1-6] show evidence for mixing 
in the D° — D° system consistent with predictions of pos- 
sible Standard Model (SM) contributions [7-11]. These 
results also constrain many new physics models [12-16]. 
An observation of CP violation (CPV) in the D° - U° 
system at the present experimental sensitivity would pro- 
vide possible evidence for physics beyond the SM [17-21]. 

One manifestation of D° — D° mixing is differing D° 
decay time distributions for decays to different CP eigen- 
states [22]. We present a measurement of charm mixing 
using the ratio of lifetimes obtained from the decays of 
neutral D mesons to CP-even and CP-mixed two-body 
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final states. We also present a search for indirect CP 
violation arising from a difference in D° and D° par- 
tial decay widths to CP-even eigenstates. Recently the 
LHCb Collaboration has reported evidence for CPV in 
the difference of the time-integrated CP asymmetries in 
D" K~K + and D° 7r-7r+ decays [23]. This mea- 
surement is primarily sensitive to direct CPV. As ex- 
plained in Appendix A, we are not sensitive to effects of 
direct CP violation at the level of the result reported by 
LHCb, and we therefore assume no direct CPV in our 
baseline model. 

We measure the effective D° lifetimes in three different 
two-body final states: K^tt^, K~K + , and 7r~7r + . We 
make no distinction between the Cabibbo-favored D° — > 
K~ir + and doubly Cabibbo-suppressed D° — > K + ir~ 
modes; in other words, we analyze and describe them 
together. Given the current experimental evidence indi- 
cating a small mixing rate, the lifetime distribution for 
all two-body final states is exponential to a good approx- 
imation. Decays in the K^n^ mode are to a CP-mixed 
final state, and are assumed to be described by the aver- 
age D° width T. The singly Cabibbo-suppressed decays 
D° (D°) to the CP-even K~K+ and 7r-7r+ final states 

are described by the partial decay rate T + (r + ), where 
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+ indicates the CP of the final state. We present in Ap- 
pendix A a discussion of the mixing formalism leading 
to the expressions that are used to extract the mixing 
parameter yep and the CPV parameter AY, 



r+ 



-r + 



yep 



AY = 



2T 



2T 



(1) 



(2) 



from the experimentally measured CP-mixed and CP- 
even lifetimes. This definition of Ay is opposite in sign 
to that in our previous measurement [2] and is now con- 
sistent with that used by the Heavy Flavor Averaging 
Group [24]. 

We use D° mesons from tagged D* + —> D°w + de- 
cays [25], as well as untagged decays where no D* + par- 
ent is found. The charge of the D*^ is used to split 
the K~K + and ir~ir + samples into those originating 
from D° and from D° mesons in order to measure the 
CP-violating parameter AY. The requirement of a D* + 
parent strongly suppresses backgrounds; hence untagged 
decays are reconstructed only in K^tt^ and K~K + be- 
cause of the relatively poor signal-to-background ratio in 
the untagged tt~tt + final state. In summary we study 
seven modes, two untagged and five tagged. 

In addition to the increased integrated luminosity of 
the new dataset compared to that used in our earlier re- 
sults [2, 3], this analysis benefits from improved particle 
identification and charged-particle track reconstruction, 
and a more inclusive and optimized event selection. We 
implement an improved background model, and we si- 
multaneously fit both the tagged and untagged datasets. 



II. EVENT RECONSTRUCTION AND 
SELECTION 



For tagged decays, we reconstruct D* + candidates by 
combining a D° candidate with a slow pion track 7r+, 
requiring them to originate from a common vertex con- 
strained to the e + e~ interaction region. We require the 
7r+ momentum to be greater than 0.1 GeV/c in the labo- 
ratory frame and less than 0.45 GeV/c in the CM frame. 
We reject a positron that fakes a 7r+ candidate by us- 
ing dE/dx information and veto any 7r+ candidate that 
may have originated from a reconstructed photon con- 
version or 7T° Dalitz decay. The distribution of the differ- 
ence Am between the reconstructed D* + and D° masses 
peaks near Am ~ 0.1455 GeV/c 2 . Backgrounds are sup- 
pressed by retaining only tagged candidates in the range 
0.1447 < Am < 0.1463 GeV/c 2 . 

To determine the proper time t and its error cr t for 
each D° candidate, we perform a combined fit to the D° 
production and decay vertices. We constrain the pro- 
duction point to be within the e + e~ interaction region, 
which we determine using Bhabha and di-muon events 
from triggers close in time to any given signal candidate 
event. We retain only candidates with a x 2 -based prob- 
ability for the fit P(x 2 ) > 0.1%, and with -2 < t < 4ps 
and a t < 0.5 ps. For tagged decays, this fit does not in- 
corporate any 7r+ information in order to ensure that the 
lifetime resolution models for tagged and untagged signal 
decays are very similar. The most probable value of tr t 
for signal events is ~ 40% of the nominal D° lifetime [27]. 

If an event contains a tagged D° decay, we exclude 
all untagged D° candidates from that event in the final 
sample. For a given final state, when multiple D° (for 
the untagged modes) or P* + (for the tagged modes) can- 
didates in an event share one or more tracks, we retain 
only the candidate with the highest P(% 2 ). The fraction 
of events with multiple D° candidates with overlapping 
daughter tracks is -C 1% for all final states. 



We use 468 fb -1 of e + e~ colliding-beam data recorded 
at, and slightly below, the T(45) resonance (e + e~ center- 
of-mass [CM] energy y/s ~ 10.6 GeV) with the BABAR de- 
tector [26] at the SLAC National Accelerator Laboratory 
PEP-II asymmetric-energy B Factory. To avoid potential 
bias, we finalize our data selection criteria, as well as the 
procedures for fitting, extracting statistical limits, and 
determining systematic uncertainties, prior to examining 
the results. 

We reconstruct charged tracks and vertices with a 5- 
layer, double-sided silicon vertex tracker (SVT) and a 
40-layer drift chamber (DCH). We select D° candidates 
by pairing oppositely charged tracks, requiring each track 
to satisfy particle identification criteria based on specific 
ionization energy loss (dE/dx) from the SVT and DCH, 
and Chercnkov angle measurements from a ring-imaging 
Cherenkov detector (DIRC). We then refit the D° daugh- 
ter tracks, requiring them to originate from a common 
vertex. To reduce contributions from D°'s produced via 
P-meson decay to a negligible level, we require each D° 
to have momentum in the CM frame pcm > 2.5 GeV/c. 



III. INVARIANT MASS FITS 

We characterize the D° invariant mass (M) distribu- 
tion for each of the seven modes with an extended un- 
binncd maximum likelihood fit to D° and D° samples. 
We allow the parameters governing the shapes of the 
probability density functions (PDFs), as well as the ex- 
pected signal and background candidate yields, to vary 
in the fits. For the tagged CP-even modes we fit the D° 
and D° samples simultaneously, sharing all parameters 
except for the expected signal and background candidate 
yields. 

We fit the tagged ir~n + invariant mass distribution in 
the fit range 1.82 < < 1.93GeV/c 2 using a sum 

of two Gaussians with independent means and widths 
for the signal PDF, along with a first-order Chebychev 
polynomial for the total background. 

The fit model for the tagged K~K + invariant mass 
distribution is similar to 7r~7r + , except that the fit range 
is 1.82 < Mkk < 1.91GcV/c 2 , and the signal PDF is 
the sum of two independent Gaussians and a modified 
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Gaussian with a power- law tail [28], which aids in better 
modeling of the lower tail of the distribution. 

The signal PDF for the untagged K~K + mode and 
for both tagged and untagged if =F 7r ± modes is a sum of 
three independent Gaussians; the background is modeled 
using a second-order Chebychev polynomial. The mass 
fit range is 1.82 < Mkk < 1-91 GeV/c 2 for the untagged 
K~K+ mode, 1.81 < M Klt < 1.92GeV/c 2 for the un- 
tagged K^ir* mode, and 1.80 < M K7! < 1.93GeV/c 2 
for the tagged if =F 7r ± mode. In these modes, we do not 
distinguish D° from D° candidates, and therefore deter- 
mine only the total signal and total background yields, in 
addition to the signal and background shape parameters. 

The reconstructed D° invariant mass distributions and 
the fit results are shown in Fig. 1, together with a plot 
of the corresponding normalized Poisson pulls [29] . 



IV. SIGNAL AND SIDEBAND REGIONS 

For the lifetime fit, we determine the regions in two- 
body invariant mass that maximize signal significance, 
minimize systematic effects due to backgrounds, and 
minimize the effect of the correlation between D° in- 
variant mass and proper time. We refer to these 
regions as the lifetime-fit mass regions. Based on 
these studies, the optimal lifetime-fit mass region is 
34MeV/c 2 wide for all tagged modes and untagged 
K^% ± events, 1.847 < M < 1.881 GcV/c 2 . Because 
of the smaller signal-to-background ratio for the un- 
tagged K~K + events, the lifetime-fit mass region for 
this mode is only 24MeV/c 2 in width, 1.852 < M < 
1.876 GeV/c 2 . For the tagged modes, a mass difference 
sideband 0.151 < Am < 0.159 GeV/c 2 is used, along with 
a low (high) invariant mass sideband, 1.819 (1.890) < 
M < 1.839 (1.910) GeV/c 2 . The low (high) mass side- 
band used for the untagged modes, 1.810 (1.899) < M < 
1.830 (1.919) GeV/c 2 , is displaced from the tagged side- 
band in order to reduce the signal component there. The 
signal purities in the lifetime-fit mass regions range from 
~ 75% for the untagged K~K + sample to ~ 99.8% for 
the tagged K^n^ events. 

We classify D° candidate decays in the lifetime-fit mass 
region as follows: D° signal decays; misreconstructed- 
charm decays, i.e., those in which the candidate- D° 
daughter tracks are decay products of a non-signal weak 
charm decay; and random combinatorial background. 
Table I gives the composition of the misreconstructed- 
charm backgrounds expected from simulated events [30] 
in each final state. 
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FIG. 1: The reconstructed two-body invariant mass distribu- 
tions for the seven modes. The vertical lines show the lifetime- 
fit mass region, defined in Sec. IV. The shaded regions are 
the background contributions. The normalized Poisson pulls 
for each fit are shown under each plot; "unt" refers to the 
untagged datasets. 



V. LIFETIME FIT 

The lifetimes arc determined from an extended un- 
binned maximum likelihood fit to t and at for candidates 
in the lifetime-fit mass region. All modes are fit simul- 
taneously using shared signal resolution function param- 



eters. The signal, misreconstructed-charm and combi- 
natorial components are described by their own set of 
PDFs, which in the tagged modes can also depend on 
the charm flavor. 

The lifetime PDF for signal is an exponential function 
convolved with a resolution function, which is the sum of 
three Gaussian functions whose widths are proportional 
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TABLE I: Expected composition (in %) of the misrecon- 
structed-charm backgrounds. Only misreconstructed-charm 
background channels that have > 1% contribution in at least 
one signal mode are listed. For the tagged modes, the yields 
are the sum of the separate D° and D° tags. 



Mode 


Ta 

7T _ 7T + 


gged Modes 
K~K+ K^-k^ 


Untagged Modes 
K~K+ 


D° Xlv 


15.4 


10.3 


29.9 


7.2 


< 2 


D° -> K--K+ 


80.8 


14.9 


57.1 


8.8 


35.8 


D° -> tt°it+K- 


1.1 


70.3 


1.7 


63.3 


6.9 


D+ -> TT + TT + K- 


< 1 


2.9 


< 1 


11.8 


< 2 


D° -> K+K~ 


< 1 


< 1 


1.3 


< 1 


3.5 


D° tt+tv- 


1.8 


< 1 


2.2 


< 1 


3.1 


D° -»• 7V+7T-7V 


< 1 


< 1 


7.0 


< 1 


17.3 


A decays 


< 1 


< 1 


< 1 


4.9 


2.6 



to at- The explicit form of the signal lifetime PDF is 
KF, L (t,a t ) = f n 'D(t,<T t ;S' T S F s 1 ,t ,T L ) (3) 
+ (!-/ti) ft2D(t,a t ;S' T S F s 2 ,t Q ,TL) 
+ (1 - ft2)T>(t,at;S^S F s 3 ,t ,T L ) , 

where fu (with i — 1,2) parameterizes the contribution 
of each individual Gaussian, Si (with i = 1,2,3) is a 
scaling factor associated with each Gaussian, and to is 
an offset of the mean of the resolution function. The 
function T>(t, a t ; s, to, r) is given by 



T>(t,a t ;s,t ,T) = 

Ca t / CXp(-i t ruc/ r ) exp 



(t-«truc+*0)' 



dt t 



(4) 



where the normalization coefficient C at is chosen such 
that 



J T>(t, at; s, to,T) dt = 1 for each <r f . 



(5) 



With this definition, the product H^[a t )-T>(t 7 a t ; s, £ , t) 
is a properly normalized two-dimensional conditional 
PDF, where £T® Ig (a t ) is a PDF characterizing the a t 
distribution, described below. To account for small dif- 
ferences in the resolution function for the different final 
states we introduce additional mode-dependent scale fac- 
tors S F , F = Kir, KK, irir. We also allow for differences 
between the resolution functions for tagged and untagged 
modes by means of scale factors S' T , T = tag (tagged) or 
unt (untagged). We fix Sk-k and S' unt to 1. 

The three lifetime parameters are tl = {t + , f + , tk^}, 
where tk-k is extracted from the tagged and untagged 
A' =F 7r ± modes, while r + and f + are extracted from the 
tagged and untagged CP-even modes. Approximately 
0.4% of the tagged CP-even samples contain correctly 



reconstructed D° candidates combined with an unre- 
lated 7r+ ; this fraction has been estimated from simulated 
events and verified in data by an earlier BABAR, analy- 
sis [1]. These candidates have the same resolution and 
lifetime behavior as those from correctly reconstructed 
D* + decays, but about half of them will be tagged as 
the wrong flavor. Therefore, the tagged CP-even D° 
proper-time distributions are modeled as the weighted 
sum of PDFs for correctly tagged and untagged candi- 
dates, characterized by the lifetime parameters r + and 
f + , respectively, and a mistag fraction /tag = 0.2%. 
The tagged CP-even D° proper-time distributions arc 
modeled in a similar fashion, where now the correctly 
tagged and mistagged PDFs arc characterized by the life- 
time parameters f + and r + , respectively. The untagged 
K~ K + proper-time distribution is modeled as a weighted 
sum of two PDFs characterized by the lifetime parame- 
ters t + and f + , respectively, and a weighting fraction 
f]ja = 0.5. These parametcrizations assume no direct 
CPV, and allow for CPV in the interference between de- 
cays with and without mixing characterized by a mode- 
independent weak phase (f>. Both / tag and / fl o are varied 
as part of the systematic error estimate for yep and Ay. 
All five tagged and two untagged signal lifetime PDFs 
arc explicitly given in Appendix B. 

The a t PDF for signal candidates is obtained directly 
from data by subtracting the sum of the background at 
distributions from that of all candidates in the lifetime- 
fit mass region. These 1-d at distributions are used to 
model the H^f{a t ) PDF discussed previously. 

We determine the t versus a t misreconstructed-charm 
signal-like PDF shape parameters and yields by fitting 
simulated events in the lifetime-fit mass region and then 
fix these parameters in the lifetime fit to data. We vary 
the lifetimes and yields as part of the study of systematic 
effects. 

The largest background in the lifetime-fit mass region 
is due to random combinations of tracks. The PDF de- 
scribing the two-dimensional combinatorial background 
in t and at in the lifetime-fit mass region is characterized 
as a weighted average of the 2-d PDFs extracted from 
the mass sideband regions. The weights for the low and 
high sidebands arc obtained from simulated events. The 
(i, at) combinatorial PDF in each sideband and for each 
mode, except for the untagged K~K + mode, is extracted 
as a 2-d histogram from the sideband samples. From 
these histograms we subtract the contribution of signal 
and misreconstructed-charm backgrounds, each of which 
is estimated from simulated events, to obtain the final 
combinatorial PDF in each sideband. For the untagged 
K~K + mode, a similar procedure is used but, instead 
of histograms, analytic signal- like PDFs are used. For 
the background PDFs the offsets and the lifetimes are 
allowed to be different for each Gaussian. The signal and 
misreconstructed-charm PDF parameters are extracted 
by fitting simulated events and then fixed, along with 
the expected candidate yields, in the fit that extracts the 
combinatorial PDFs in each sideband. 



S 



For the untagged K~K + mode both the expected 
signal and combinatorial yields are free parameters in 
the lifetime fit. The expected combinatorial background 
yields in the other modes are determined by integrat- 
ing the total background PDF extracted from the mass 
fit in the lifetime-fit mass region, and then subtracting 
the expected misreconstructed-charm background yields, 
which are determined from samples of simulated events. 
A small bias on these fit yields is observed in fits to simu- 
lated events. To correct for this, we scale the data yields 
based on the simulated-event fits and vary the mode- 
dependent scale factors as a systematic uncertainty. Ta- 
ble II gives the event class yields plus uncertainties ob- 
tained from the lifetime fit and indicates the yields that 
are fixed. 

TABLE II: Signal and background yields in the lifetime-fit 
mass region. Yields with uncertainties are those obtained 
directly from the lifetime fit to data. For the tagged modes, 
the yields are the sum of the separate D° and D° tags. 

Tagged Modes Untagged Modes 

7T~7r+ K~K+ Jf T 7T ± K~K+ A' T 7T ± 

Signal 65 430 136 870 1487 000 496 200 5 825 300 

±260 ±370 ±1200 ±1200 ±2600 
Comb. Bkgd. 3760 653 2849 165 000 1044 552 

±1000 

Charm Bkgd. 97 309 642 5477 4645 



The simultaneous fit to all events in the lifetime-fit 
mass region has 20 floating parameters: the seven signal 
yields and three signal lifetimes; the yield of untagged 
K~K + combinatorial candidates; the offset to; the pa- 
rameters fn and ft2 characterizing the weight of each 
Gaussian in the signal resolution mode; and the proper- 
time error scaling parameters s±, S2, S3, Skk, SW, and 
S[ ag . After extracting the three signal lifetimes, using 
their reciprocals in the computation of yep and Ay as 
defined in Eqs. (1) and (2), respectively, we find 

y C p = [0.72±0.18(stat)]%, 
AF = [0.09±0.26(stat)]%. 

The statistical errors are computed using the covariance 
matrix returned by the fit. The lifetime-fit mass region 
proper-time distributions and projections of the lifetime 
fit for the seven different decay modes are shown in Fig. 2. 



VI. CROSS CHECKS AND SYSTEM ATICS 
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FIG. 2: Proper-time t distribution for each decay mode with 
the fit results overlaid. The combinatorial distribution (in- 
dicated as 'Comb.' in light gray) is stacked on top of the 
misreconstructed-charm distribution (indicated as 'Charm' in 
dark gray). The normalized Poisson pulls for each fit are 
shown under each plot; "unt" refers to the untagged datasets. 
The bottom right plot shows the individual lifetimes (with 
statistical uncertainties only); the gray band indicates the 
PDG D° lifetime ±lcr [27]. 



We have performed numerous cross checks to search 
for potential problems, in addition to quantitative studies 
that yield the systematic uncertainties given in Table III, 
discussed below. Initially we tested the fit model by gen- 
crating large ensembles of datasets randomly drawn from 



the underlying total PDF, and observed no biases in the 
yep and Ay results obtained. In addition, we have fit 
an ensemble of four simulated datasets, each equivalent 
in luminosity to the data, and found no evidence of bias 
m Ucp or Ay. 
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In fitting the data, we find that the tagged and un- 
tagged extracted lifetimes for K~ K + , and separately 
for K^n^. arc compatible within the statistical uncer- 
tainties. We performed a simultaneous fit to the tagged 
channels, and a separate simultaneous fit to the untagged 
channels, and find the lifetimes to be compatible within 
the statistical uncertainties. We repeated the fit allow- 
ing the K~K + and ir~ir + final states to have separate 
t + and f + lifetimes, and observed no statistically signif- 
icant difference between the K~K + and -k~it + results. 
We estimated the effects of the SVT misalignment to be 
negligible. 

We varied the lifetime-fit mass region width by ±4 and 
±2 MeV/c 2 . We adopt as the systematic uncertainty half 
the RMS of the differences | A[y C p] | and | A[AY] | from the 
nominal fit central values. We also shifted the position 
of each mass region by centering each of them at the 
most probable value for the signal PDF obtained in the 
invariant mass fits. These systematic uncertainties are 
given in the first two lines of Table III. 

For the untagged K~ K + mode, the combinatorial 
yield is a parameter determined in the lifetime fit. How- 
ever, it is also needed to determine the signal a t PDF. We 
first use the total background yield determined from the 
mass fit to extract a signal a t PDF, which is employed 
in an initial simultaneous lifetime fit. The combinatorial 
yield from this fit is used to construct an improved at 
signal PDF and a second fit is performed (the nominal 
fit). We estimate the systematic error on yep and AY 
associated with the determination of the signal at PDF 
for the untagged K~K + mode to be the difference in the 
values obtained from an additional iteration of the fit and 
the nominal fit. 

We vary the nominal mistag rate of 0.2% by ±0.04%, 
a 20% relative variation, and find no significant change 
in the nominal fit values. Instead of assuming equal frac- 
tions of D° and D° in the untagged K~K + mode, we 
adopt the latest CDF result for direct CPV [32] , and find 
negligible change in yep and AY. 

We rely on simulated events to determine both the 
PDF shapes and yields for the misreconstructed-charm 
backgrounds. To account for the model dependence, we 
vary the effective lifetime of these events by ±5%, ex- 
cept for the tagged 7r~7r + mode where the variation is 
±15% due to the small number of simulated events that 
pass the selection criteria for this mode. We also vary 
the expected misreconstructed-charm yields by ±10% in 
the tagged channels, and ±5% in the untagged channels. 
Each variation is simultaneously applied to all modes. 
These are > 2a variations relative to the statistical un- 
certainties of the simulated datasets. 

We vary the yields, weighting parameters, and fit- 
ting strategy used to obtain the 2-d lifetime PDF for 
combinatorial-background events in the lifetime-fit mass 
region from the mass sidebands. The yields for the tagged 
combinatorial-background events are varied by ~ 5% in 
the 7r _ 7r + mode, 15% in tagged K~K + , and 20% in 
A' =F 7r ± . The untagged AT =F 7r ± combinatorial-background 



TABLE III: The yep and AY systematic uncertainties. The 
total is the sum-in-quadrature of the entries in each column. 



Jj 11 V dlleLLrlUIl 


\^[ycp\\ \/o) 


1 ATA VII (°&\ 


mass window width. 


U.UO ( 


u.uzz 


mass window position 


n nn^ 
u.uuo 


n nm 

U.UUl 


untagged KK signal o± PDF 


n D99 
u.uzz 


n nnn 
u.uuu 


mistag fraction 


n nnn 
u.uuu 


n nnn 
u.uuu 


untagged KK D fraction 


n nm 


n nnn 
u.uuu 


charm bkgd. lifetimes 


n pi/1 o 
U.U4Z 


n nm 
U.UUl 


charm bkgd. yields 


0.016 


0.000 


comb, yields 


0.043 


0.002 


comb, sideband weights 


0.004 


0.001 


comb. PDF shape 


0.066 


0.000 


at selection 


0.052 


0.053 


candidate selection 


0.028 


0.011 


Total 


0.124 


0.058 



yield is varied using the value extracted from an alter- 
native lifetime-fit model in which the yield is allowed to 
vary. 

The weights given to the low- and high-mass sidebands 
in the data in order to derive the combinatorial PDF in 
the lifetime-fit mass region in data are extracted from 
simulated events. They are varied by plus and minus 
the statistical uncertainty derived from splitting the sim- 
ulated dataset, which is equivalent to several times the 
nominal integrated luminosity, into datasets that numer- 
ically match the nominal luminosity. 

We also apply the variations described above for the 
misreconstructed-charm background to vary the yields 
and shape of the PDF that describe the residual sig- 
nal events in the sidebands. This is also done for the 
misreconstructed-charm PDF used in the sideband fits 
from which the 2-d combinatorial PDF is extracted. This 
yields the combinatorial PDF shape variation, which is 
then used in the nominal fit, to obtain the variation re- 
ported in Table III. 

Finally, we vary the at criteria by ±0.1 ps from the 
nominal at < 0.5 ps, and take as the systematic uncer- 
tainty the RMS of the deviations from the nominal fit 
central value divided by \[2. We also consider two vari- 
ations in how multiple candidates sharing one or more 
daughter tracks are treated. In the first variation, we re- 
tain all multiple candidates, if each candidate passes all 
the other selection criteria. In the second variation, we 
reject all multiple candidates sharing one or more daugh- 
ter tracks. We fit these datasets using the nominal fit 
model, and assign the largest observed deviation from 
the nominal yep and AY" central values as the system- 
atic uncertainty in Table III. The total yep and AY 
systematic uncertainties are calculated by summing the 
contributions from all sources in quadrature, and are re- 
ported in the last row of Table III. 
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VII. CONCLUSIONS 
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Appendix A: Mixing Formalism and Considerations on the Role of Direct CP Violation 



The time evolution of the flavor eigenstates D° and D° is governed by the Schrodinger equation: 




The mass eigenstates D\ and D 2 are obtained from the diagonalization of the effective Hamiltonian H c s — M — |I\ 
Under the hypothesis of CPT conservation the two mass eigenstates can be written in terms of the flavor eigenstates 
as 



D°(t) 



|Z>!) = p\D°}+q\D ) , 
\D 2 ) = p\D°)-q\D°) , 



(A2) 



where 



2 



M12 - |r 12 



and |p| 2 + |g| 2 = 1. (A3) 



We choose the positive root for q/p; choosing the negative one just means exchanging D\ with D 2 . If CP\D°) = +\D°), 
in case of no CSV, D\ is the CP-even state and D 2 the CP-odd state. 

It is traditional to quantify the size of D° — D° mixing in terms of the parameters x = Am/T and y = Ar/2r, where 
Am = mi — 7712 ( Ar = Ti — T 2 ) is the difference in mass (width) of the states defined in Eq. A2 and V = (Ti + T 2 )/2 
is the average width. If either x or y is non-zero, mixing will occur. While most Standard Model expectations for the 
size of both are < 10~ 3 [10, 33], values as high as 10 -2 or even higher are predicted by certain models [13, 15]. 

CP violation can manifest in D° decays in three ways: 

• in decay, when \Af/Af \ 7^ 1, 

• in mixing, when r m = \q/p\ =/= 1, 

• in the interference between decays with and without mixing, when the weak phase 0/ofA/ = ^-jjis different 
from zero, 

where A f (A f ) is the amplitude for D° (D°) decaying into a final state /, A f = (f\H D \D°) (A f = (f\H D \D )). 
The presence of mixing alters the exponential distribution for the D° decay into a final state /. In particular we 
have 

T(D°(t) -)■ /) = §\A f \ 2 e- n [(l + |A / | 2 )cosh y ri+(l- |A/| 2 ) cosxTi (A4) 

-23?(A / )sinhyrt + 23(A / )sina;R] , 

r(S°(t) -> /) = \\A } \ 2 e- Tt [(l + |A7 1 | 2 )cosh y rt+(l-|A7 1 | 2 )cosxn (A5) 

-2SR(A7 1 )sinh y n + 23(A7 1 )sina;r< . 

In this analysis we are interested in CP-even final states (/ = h + h~ ,h = K, n). If we neglect second-order terms 
in xTt and yTt, the decay time distributions can be treated as exponentials with effective widths [34]: 

r(D°(()^/) « e< f with r+ h = r [1 + y $t(\ hh ) - x %(\ hh )} , (A6) 

T(D°(t) /) cx e^- 4 with Tt lh =r[l + yM(\-l)-xZ(\ M l)}. (A7) 

To better understand the effects of CP violation we introduce two more parameters, one describing CfV in decay 
(Ap) and one in mixing (Am)'- 



\A f /A f f-\A f /A f f 
Ad \A f /A f \* + \A f -/A f -\*> (A8) 



Am - 7 , - a - (A9) 
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Since / = h + h then / = /. Noting that there is no strong phase in A/ since the final state is its own CP-conjugate, 
we can express A n „ in terms of A^ 1 , Am and the CP- violating phase <$>hh'- 



1 - i'n 1 + A 



M 



1 + A% I -A M 



1/4 



(A10) 



Expanding Eqs. A6 and A7, and retaining only terms up to first order in A^y and Am, we obtain 



r+ 

1 hh 



1 



1 + (y cos (fthh - x sin (f>hh) + T^i^M - Ad ) (j/ cos ^ _ x sin ) 



AmAjj (y cos - a; sin <^h) 



(All) 



r + 



1 + (y cos + x sin - -(4 M 



--A M Ap l (y cos 0hh + x sin 0^,) 



(y cos ^ + xsin(j) hh ) 



(A12) 



Combining the widths defined above we obtain the two observables yep and Ay, which, in general, depend on the 
final state because of the CPV parameters A 1 ^ and cf>hh- 



Vcp ~ 



r 



ill i 



r 



hh 



2T 



\yhh _ hh ran 

2r 



Other experiments characterize the CP-violating observable as Ar, 



r+ 



hh 



r hh r hh 



The relationship between A-p, AY and yep is 



Ay = (l+ycp) A r . 



(A13) 
(A14) 

(A15) 
(A16) 



These quantities are directly related to the fundamental parameters that govern mixing and CPV in the charm 
sector: 
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A M A h £ y cos <j) hh , 



(A17) 
(A18) 



Both yep and Ay are zero if there is no D° — D mixing. Otherwise, a non-zero value of yep implies mixing and a 
non-zero value of Ay implies CPV . 

In the charm sector, because the CKM elements involved belong to the Cabibbo submatrix, we can assume that 
the weak phase <j>hh does not depend on the final state: 4>hh — 4> [35]. As stated earlier if direct CPV has a significant 
effect, then the values of yep and Ay depend on the final state. In this analysis we assume that the effect of direct 
CPV is negligible in the decays to CP eigenstates; i.e., we assume = r+ ff (and T^ K = T^ w ). In Eqs. All 

and A12 this means neglecting the linear terms in A h ^ . Assuming that Ajj and y are both 0(1%) and 4>hh = 0, the 
neglected term is 0(1O -4 ), beyond any current experimental sensitivity. 

Under the above assumptions, Eqs. All and A12 become 
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r + 



1 + (y cos 4> — x sin < 
1 + (y cos 4> + x sin < 



(y cos <j) — x sin < 
(y cos <j) + x sin < 



(A19) 
(A20) 
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Inserting Eqs. A19 and A20 into Eqs. 1 and 2 yields 

Vcp = ycos0 - ^-xsincj) , (A21) 
Ay = -xsm(j>+ ^-ycoscfi. (A22) 

From the experimental point of view, we measure three lifetimes instead of the partial widths: 

• r + for the D° — s- K~K + , ir~n + decays, 

• f + for the D° —> K~K + . ir~ir + decays, 

• for the D° (and D°) K^ir decays (the Cabibbo favored K~tt + and the doubly Cabibbo suppressed 
K + ir~ decays are collected in the same sample), 

and use their inverse to compute yep and AY. 

The measured observablcs constrain the parameters that govern mixing and indirect CPV in the charm sector. 

Appendix B: Signal Lifetime PDFs 

The explicit form of the signal lifetime PDFs based on the prototype PDFs presented in the main text are given 
below: 
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where f^„ = 0.2%, /po = 0.5 and Sk-k = S' unt = 1 are fixed in the nominal fit. 
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